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Abstract 

In this article we define Kober fractional integral operators in the multivariable case. First we consider 
one sequence of independent random variables and an arbitrary function, which can act as the joint density 
of another sequence of random variables. Then we define a concept, analogous to the concept of Kober 
operators in the scalar variable case. This extension is achieved by using statistical techniques and the 
representation gives an interpretation in terms of a joint statistical density. Then we look at two sets of 
random variables where between the sets they are independently distributed but within each set they are 
dependent. Again extensions of Kober fractional integral operator are considered. Several such statistical 
interpretations are given for Kober operators in the multivariable case. 

1. Introduction 

When going from a one-variable function to many-variable function there is no unique one to one 
correspondence. Many types of multivariable functions can be considered when one has the preselected 
one-variable function. Hence there is nothing called the multivariable analogue of a univariable operator. 
Hence we construct one multivariable operator here which is analogous to a one variable Kober fractional 
integral operator of the second kind. Other such analogues can be defined. 

Definition 1.1. Kober fractional integral operator of the second kind in the multivariable case This will 
be defined as the following fractional integral and denoted as follows: 

3=1 V 31 



k 

X 



{f[f (« i -« J -) a, "S C ' aj }f(vu...,v k )d Vl A...Adv k . (1.1) 



This definition is parallel to the one in the one variable case. We will now look at various connections 
to different problems. First we will establish a number of results in connection with statistical distribution 
theory We will show that (1.1) can be treated as a constant multiple of a joint density of a number of 
random variables u\, u k appearing in different contexts. 

1.1. Kober fractional integral operators of the second kind in multivariable case as statistical 
densities 
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Let xi, X2, Xk be independently distributed typc-1 beta random variables with parameters (Q + 
1, dj),j = 1, k, Q > — 1, a.j > 0, j — 1, k. Usually the parameters in a statistical density are real but 
the following integrals also exist for complex parameters and in that case the conditions will be 9?(£) > — 1 
and $l(aj) > 0, j = 1, k. That is, the density of Xj is of the form 

= Kt!! '^ 1 -^" 1 ' • • 1 (1-2) 

for Qfj > 0, Q > —1 and fj(xj) = elsewhere, j = 1, k so that the joint density of x\, x k is the product 
fi(xi)...fk(xk)- Let Vk be another sequence of real scalar positive random variables having an arbitrary 

joint density f(v\, Vk). Let the two sets (xi, Xk) and (v\, Wfe) be independently distributed. Consider 
the transformation Uj = XjVj or Xj = ^f- and the Jacobian of the transformation is 

{f[ Adz,} A {f[ Ad^} = {fliv^ufl Ad Uj } A {f[ Ad Vj } (1.3) 
i=i j'=i j=i j=i j=i 

Then the joint density of Hi, Mfc, denoted by Ufe), is given by 



x /(wi,..., w fe )di; 1 A... A dw fe 

x /(ui, ...,w fe )dwi A ... Adwfe. (1.4) 

Therefore one can write 

Theorem 1.1. Let Xj,Uj,Vj,j — l,...,k be as defined above where Xi,...,Xk are independently type-1 beta 
distributed with parameters ((j + l,aj),j = l,...,k, vi,...,Vk having a joint arbitrary density f(v\,...,Vk) 
with (x\,...,Xk) and (vi,...,Vk) are independently distributed. If the joint density of u\, ...,Uk is denoted as 
g(ui, -,u k ) then 

(fl rfo + + 1) }9{UU Uk) = -, u k ) (1.5) 

for M(( j )>-l,$l(a j )>0,j = l,...,k. 

In this case we had x\,...,Xk mutually independently distributed and thus there were a total of k + 1 
densities involved, the k of Xi, ...,Xk and the one of (v\, ...,Vk)- Let us see what happens if X\, ...,Xk are not 
independently distributed but they have a joint density fi(xi, ...,£&) and (v\, ...,Vk) having a joint density 
f2(v i, Vk). Then we can show that if /i can be eventually reduced to independent type-1 beta form, 
still we can consider Kober operators of the second kind in the multivariable case as constant multiples of 
statistical densities. 

Let (xi, Xk) have a joint type-1 Dirichlet density with parameters (ai + 1, a k + 1; afc+i), $l(aj) > 
-1,3 = l,...,k,3t(a k +i) > 0, that is, 



h(xi,..,xk) = r i ai + - + Qfc+1 + fc) xT...xT 
{n- =1 r(a, + i)}rK +1 ) 

x (1 -xi - ...-xk) ak+1 ~\0< xj < 1,3 = l,...,fc,0 < xi +... +x k < 1 (1.6) 
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and fi(xi, Xk) — elsewhere. Let us consider the transformations x\ = yi,x 2 = 2/2(1 — yi),... Xk = 
(1 - yi)...(l - y k -i) or 



x 3 = vA 1 - yOC 1 - y2)-(i - yj-i),j = 

or 

^' l-^-'-^ -^ 1 -*- (L7) 

Under this transformation the Jacobian is (1 — j/i) fe_1 ...(l — j/fc-i). It is easy to show that under this 
transformation yi,...,yk will be independently distributed as type-1 beta variables with the parameters 
(ctj + 1, (3j) with /3j = ctj+i + aj+2 + ■■■ + ctk + (k — j) + ctk+i or yj has the density 

and fj(yj) = elsewhere, a, > > 0,j = 1, k. In the light of these observations, let us consider two 

sets of positive random variables (xi, ■■■,Xk) and (vi,...,Vk) where the two sets are independently distributed 
with (xi, Xk) having a type-1 Dirichlet distribution. Let Uj = yjVj = Vj( 1 _ Xl J C3 — . ),j = 1, k. Then 
following through the same procedure as above we have the following theorem. 

Theorem 1.2. Let (xi,...,Xk) and (vi,...,Vk) be two sets of real scalar positive random variables where 
between sets they are independently distributed. Let (vi,...,Vk) have an arbitrary joint density f(v\,...,Vk) 
and let (x\, ...,Xk) have a type-1 Dirichlet density with the parameters (ct\ + 1, ...,ak + l;afe + i) or with the 
density 

fi{xi, ...,x k ) = C X? ...x^(l - Xl - ... - x^- 1 (1.8) 
< xj < 1, < x\ + ... + Xk < l,j = 1, k 

and fi(xi, ...,Xk) = elsewhere, where C is the normalizing constant. Letuj = Vj( 1 _ xi J Cj _ x ._ 1 ),j = 1, k. 
If the joint density of u\, ...,Uk is again denoted by g(ui, ...,Uk) then 

in ffcj + 1 + 1) >g^> ■•'"*)= <tl j :t- k f^u ..,u k ) (1.9) 

where Pj = a j+1 +a j+2 + ■■■ + a k + (k-j) + a k +i,j = l,...,k, $l{aj) > -l,j = 1, k, $t(a k+ i) > 0,SR(/3j) > 
0,j = l,...,k. 

The above structure indicates that we can consider any multivariable density fi(xi, ■■■,Xk) for a set of 
real scalar positive random variables {x\ , Xk) and if we can find a suitable transformation to bring the joint 
density of the new variables as products of type-1 beta densities then the Kober fractional integral operator 
of the second kind for the multivariable case can be written in terms of a statistical density as shown above. 
There are many densities where a transformation can bring f\{x\, .., Xk) to product of type-1 beta densities. 
There are several generalizations of type-1 and type-2 Dirichlet models where suitable transformations exist 
which can bring a set of mutually independently distributed type-1 beta random variables. We will list one 
more example of this type before quitting this section. 

Let us consider a generalized type-1 Dirichlet model of the following type. Several types of generalizations 
of the following category are available. 



h( Xl , x k ) - d x^ (1 - Xl )^x^(l - Xl - x 2 f\.. 

x xl k {l - Xl - ... -xkf k ~\0 < x 1 + ... + Xj <l,j = l,...,fc (1.10) 
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and fi(xi, Xk) = elsewhere, where C\ is a normalizing constant. Let us consider the same transformation 
as in (1.7). Then we can show that yi, ...,yk will be mutually independently distributed as type-1 beta random 
variables with the parameters (ctj + — 1, fc where 

7j = ctj+i + a j+2 + ... + a k -i + fa + Pj+i + ••• + Pk + (k - j) (1.11) 

for j = 1, fc with 5i(o!j) > — 1, j = 1, fc and 3?( / y J ) > 0, j = 1, fc. 

Theorem 1.3. Let Xi, ...,Xk have a joint density of the form in (1.10). Let v\, ..,Vh be another set of real 
scalar positive random variables having an arbitrary density f(vi,...,Vk). Between sets let (xi,...,Xk) and 
(vi,...,Vk) be independently distributed. Consider the transformation as in (1.7) where 

Let the joint density of u\, ...,Uk be again denoted by g(u\, Uk). Then 

r( aj + 1) 



*/(»!,....»*) (i.i2) 

where jj = aj + i+a j+ 2+...+ak+Pj+Pj+i+.-+Pk+(k-j),j = 1, fc for 3? (ay) > -1, Sft(7j) > 0, j = 1, fc. 

1.2. A Pathway Generalization of Kober Operator of the Second Kind in the Multivariable 

Case 

Let xi,...,Xk be independently distributed with Xj having a pathway density given by 

fi( Xj ) = c jp xf [1 - a,(l - qfixjX&i (1.13) 
for 1 — Oj (l — Qj)xj > 0, dj > 0, < 1, 77^ > 0, Cj > —1 and fj(xj) = elsewhere, where 

Mi -^o+^fe + i + ^ + i) 

rto + iM^ + i) (L14) 

Let be real scalar positive random variables with a joint density f(vi,...,Vk). Let (xi,..,Xk) and 

(til, Vk) be statistically independently distributed. Let Uj — XjVj, xj = ^f, j = 1, fc. Then the Jacobian 
of the transformation is {v\...Vk)~ 1 - Let the joint density of u\, Uk be denoted by g(ui, Uk). Then from 
the standard technique of transformation of variables the density g is given by 



/>OC 

g(u-L,...,u k ) = !l{ <%,.,»; / 



-0-(t%+D 

• '/./)" 
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x [vj - a(l - qj)uj] 1 ^ }/(«!,..., w fe )dw 1 A ... A dw fe . (1.15) 
Hence we may define a pathway extension of Kober operator. 

Definition 1.2. 4 Pathway Kober Fractional Integral Operator of the Second Kind for the Multivariable 
Case It will be defined and denoted as follows: 

(C- Vj I 1) 7=1 k k u'' f n 3 

j=l Vl-gj " r 1 / Jv j >a j (l-q j )u j 

x /(wi, ...,v k )dv! A ... A dufe. (1.16) 



for a,j > 0,qj < l,T]j > 0,K(Cj) > -1. 

Theorem 1.4. Let Xi, Xk, v\,...,Vk, Uj,j — l,...,k and g{u\, ...,Uk) be as defined in (1.15). Let the 
pathway extended Kober operator be as defined in (1.16). Then 



(1.17) 



When any particular g r — ► 1_ then we can see the corresponding factor going to the exponential form. 

(i-g r )Cr+ir(c r + T ^. + i) ! 



lim 

9r-»-i- 



IW- + 1) 



-(^-[l-a r (l -<?,)(- 



l-9r 



= (!t)Crl e -r^ (^) j0 <« r <OO. 



(1.18) 



Thus, individual ^'s can go to 1 and the corresponding factor will go to exponential form or the correspond- 
ingly we get a gamma density structure for that factor. 

1.3. Mellin Transform in the Multivariable case for Kober Operators of the Second Kind 

The Mellin transform in the multivariable case is defined as 

POO poo 

M{f(x u ...,x k );s u ...,Sk}= ... xl 1 - 1 ..^- 1 

Jo Jo 

x f(xi, ...,Xk)dxi A ... A dx fe (1.19) 

whenever it exists, where si, Sk in general are complex parameters. Hence for the Kober operator of the 
second kind we have 

pOO poo 

•W/C; ; :l.,! ./•>": «*} I ...J o «• 1 



j=l ^( a o) JVj>U 3 

x f(v u ...,v k )dV dU 
where, for example, dU — dui A ... A duk 



• I II II ^— / (^""iP V " J \ 



/ -jf M ''^11'^ " 'III./, »?' Cj '(''./ -»i)* >J l <^:i}W- 



Take out Uj, put yj = ^ then the integral over will go to 



rK + Ci-.si)' 



Hence the required Mellin transform is 



r(Q + gj) 
M[ r ("j + + 



n 



/*(S1, ...,Sfe) 



where /* is the Mellin transform of /. Then we have the following theorem. 
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Theorem 1.5. For the Kober fractional integral operator of the second kind in the multivariable case as 
defined in (1.1) the Mellin transform, with Mellin parameters si, ...,Sfc is given by 



^{4?j=i^'' ,fc /(«i,-.-,«fc);*i,-»,«fc} 



LJ-r^ + o + s,) 



n 



/*(si,...,Sfc) 



(1.20) 



for St( aj ) > 0, fft(d + sj) > 0,j = 1, fc. 

2. Kober Fractional Integral Operator of the First Kind for Multivariable Case 

In the multivariable case we will start with the following definition and notation. 
Definition 2.1. Kober Fractional Integral Operator of the First Kind in the Multivariable Case 



^nnV/ 

7=1 ^ 3> J Vj=Q 



(uj-vjr^vf} 



x f(v u ...,v k )dV. 



(2.1) 



First, we will derive this operator as a constant multiple of a statistical density. To this end, let x\, Xk 
be independently distributed type-1 beta random variables with the parameters (Q,etj),j — 1, k or with 
the density 



^) = f^f^ !B ?" 1 ( 1 -^ a, - 1 ' <^< 1 ' 



(2.2) 



for ctj > 0, Q > or 5J(ay) > 0, $t(Cj) > when the parameters are in the complex domain. Let (i>i, Vk) 
be real scalar positive random variables having a joint density f(v\,...,Vk)- Let Uj = ^f-,j = l,...,k. The 
Jacobian is given by 

k 

dX AdV =([[(-^)}dU AdV (2.3) 



where the earlier simplified notation is used. The joint density of U\,...,Uk, following through the earlier 
steps, denoted again by g(ui, Uk), is given by 



x f( Vl ,...,v k )dV. 



(2.4) 



Hence we have the following theorem. 

Theorem 2.1. Let Xi,...,Xk be independently distributed type-1 beta random variables and let Vi,...,Vk, 
u\, ... , Uk be as defined in (2.2) and (2.3). Let the joint density o/mi, Uk be denoted by g(u\, Uk). Then 



(2.5) 



We can have pathway extension to Kober operator of the first kind, parallel to the results for the case 
of second kind. Other properties follow parallel to those for the case of the operator of the second kind. We 
will evaluate the multivariable Mellin transform following through steps parallel to those in the case of the 
second kind and hence we give the result here as a theorem. 
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Theorem 2.2. The Mellin transform, with Mellin parameters si, Sk, for Kober fractional integral oper- 
ator of the first kind in the multivariable case is given by the following: 

k 

M{K^fjl : j :i- k f(u u u k ); Sl , s k } = {H _JXl±<L_fL_ }/ * (si; Sfe) (2 . 6) 

for 3i(aj) > 0, SR(Cj) > 0, Sft(s) < 1 + 3ft(Cj ) > J = 1) ■••> fc w^ere /* is tfte MeHin transform of f. 

We can obtain theorem parallel to the ones in section 1. Some of these will be stated here without 
proofs. The derivations are parallel to those in section 1 and hence omitted. 

Theorem 2.3. Let x\, x k be independently distributed as the pathway model in (1-13) with Q replaced by 
Q — 1, j = 1, k. Let (vi,....,Vk) have a joint arbitrary density f(v±, ...,v k )- Let (xi,...,Xk) and (vi, v k ) 
be independently distributed. Let Uj — or xj = ^- ,j = 1, k. Let the joint density of U\, ...,Uk be again 
denoted by g(ui, ...,u k ). Then 

n ■ k _ O ~ ( 1-q ■ + "0 u j 

c;.^-> k «*) = in n_a_, n / - - ^ ^ > 

x f{v 1 ,...,v k )dv 1 A ... Adv k (2.8) 

/or a, > 0, qj < l,r)j > 0,R(Cj) > 0,j = l,...,fc. 

From here we can have a definition for a pathway extension of Kober operator of the first kind in the 
multivariable case. 

Definition 2.2. Kober fractional integral operator of the first kind in the multivariable case Let the vari- 
ables and parameters be as defined in Theorem 2.3. Then the pathway extended Kober operator of the first 
kind is defined and denoted as follows: 



l u J .„ J .^. j -i a- -. u k) = {[[ jT, a , u / \ u i - a A l - 1-93 V 1 

x f(vi,...,Vk)dvi A ... A dv k 

for aj > 0,gj < l,rjj > 0, Ji(Cj) > 0,j = l,...,fc. 

We can list several theorems when x\, ...,x k are not independently distributed. Two such cases will be 
listed here without proofs. The proofs will be parallel to those in section 1 and hence omitted. 

Theorem 2.4. Let xi,...,x k have a type-1 Dirichlet density as in (1.6) with aj replaced by aj — 1 for 
j = l,...,k and let (vi,...,v k ) have an arbitrary joint density f(vi,... 7 v k ) where (xi,...,x k ) and (vi,...,v k ) 
are independently distributed. Let yj = 1 _ Xi _ Xj _ x . i and let yj = ^f-,j = l,...,k. Let the joint density of 
U\, ...,u k be denoted by g(ui, u k ). Let (3j = ctj + i + aj +2 + ... + a k . Then 

^^'•"•v^i, ••-,«*) = {fi f^fk) }g(Ml, - ,Mfc) - (2 - 9) 



The next theorem is parallel to theorem 1.3 for the operator of the second kind. 

Theorem 2.5. Let Xi,...,x k have a joint density as in (1.10) with aj replaced by aj — l,j = l,...,k. 
Let yj = 1 _ Xl J Ll — - as defined in (1.7). Let (vi,...,v k ) have an arbitrary joint density f(vi,...,v k ). 
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Let (xi, Xk) and {v\,...,Vk) be independently distributed. Let Uj = ^ or yj = = l,...,k. Let 

7j = a J+ i + a i+2 + ... + a k + /3j + + ... + (3 k . Then 

4gS = {fl r( ^ ( " j) 7j) }g(t«i,...,« fc ) (2-10) 

/or K(aj-) > 0,5R(7j) > 0, j = l,...,k. 
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